:arXiv:nlin/0009020vl [nlin.CD] 8 Sep 2000: 



CD 

5' 
P 

cc 
-: 

r+ 
5' 



Asymptotically exact spectral estimates 
for left triangular matrices 

Michael Blank* 

Russian Ac. of Sci., Inst, for Information Transmission Problems, 
B.Karetnij 19, 101447, Moscow, Russia, blank@iitp.ru 

September 7, 2000 



Abstract 

For a family of n * n left triangular matrices with binary entries 
we derive asymptotically exact (as n — ► oo) representation for the 
complete eigenvalues-eigenvectors problem. In particular we show that 
the dependence of all eigenvalues on n is asymptotically linear for large 
n. A similar result is obtained for more general (with specially scaled 
entries) left triangular matrices as well. As an application we study 
ergodic properties of a family of chaotic maps. 
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1 Binary left triangular matrices 

By an n x n left triangular matrix with binary entries we mean the matrix 
A n = (dij) whose entries up to (and including) the secondary diagonal 
are ones, while all others are zeros, i.e. aij = 1 if and only if 1 < i < 
n — j + 1. Despite its very classical appearance spectral properties of such 
matrices were not known, probably due to the fact that for a finite matrix 

*This research has been partially supported by INTAS, CRDF and RFBR grants. 



size n there is no reasonable representation for the spectrum. Denote by 
A} , Ajj > ■ ■ ■ j ^n the eigenvalues of the matrix A n ordered by their modules 
(in the decreasing order) and by {e k n }J? =1 the corresponding orthogonal 
system of eigenvectors (which exists due to symmetry of the matrix). As 
usual by (e), we denote the i-th entry of the vector e G R ra and by ||e|| - its 
L 2 -norm. 

We start from numerical results. It turns out that numerically with 
very high accuracy for each fixed k £ {1, . . . ,n} the dependence of the k- 

(n) 

th eigenvalue X k is linear on n. From this result one might expect that 
the determinant det(^4 n ) grows as n n for large n, however expanding the 
matrix with respect to the last line one shows that det(.A n ) = (— l) n+1 . 
So there is no surprise that another numerical test shows that for a fixed 

(n) 

value n the modules of eigenvalues |AJL | decay hyperbolically on k. The 
following statement gives an asymptotic representation for the spectrum, 
which confirms these numerical findings. 

Theorem 1.1 Let A n be the left triangular matrix. Then for each fixed 
k£{l,...,n} 

(«!•>), = m ( <*- ''f - 1) ,) + 0(1/ ,o. 



The proof is based on the following statements. First notice that 

n- 

r- 



(A n v)i = YJj={ +1 v j f° r an Y vector v G R n . 



Lemma 1.1 For any C 1 function f : [0, 1] — ► [— 1, 1] with the derivative not 
exceeding 1 by absolute value we have 



— i-hl Til — iH - ! i i i — i 

~1 ~i n Jo 



n—i+l n—i+1 

E m - 1)/«) = n e 

Proof. Direct computation. 



Lemma 1.2 The operator L : f(x) — ► Jq~ x f(s) ds considered as an opera- 
tor acting in the space L 2 has the complete system of orthogonal eigenfunc- 
tions 

E k (x) := cos({k - l/2)irx) 



with eigenvalues 






Vk ■= ( 


fork = 1,2,.... 




Proof. Indeed, 
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cos((/c - 1/2)ttx). 
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Lemma 1.3 Lei i/ie equality Av = /it>+£ &e satisfied for a symmetric matrix 
A, two vectors v,£ € R n with ||£|| < e[|u|| and a scalar fj,. Then the closest 
to fj, eigenvalue A of the matrix A satisfies the inequality jA — fi\ < e and 
if its multiplicity is equal to 1, then the corresponding eigenvector e is such 
that \\e — v\\ < 0(e) ■ \\v\\. 

Proof. Probably this statement is not quite new (see for example close 
statements in |J), however since we were not able to find exact refer- 
ences to needed estimates we give a sketch of the proof. Notice that (A — 
£,v* /\\v\\ 2 )v = [j,v, i.e. the vector v is the eigenvector of the perturbed matrix 
A^:= A- £v* /\\v\\ 2 . Indeed, 

Afj,v = [iv + £, — £,v*v/\\v\\ 2 = [iv + £ — £ = fiv. 

On the other hand, 

[|^7IMI 2 ||<^<£. 

Thus the matrix A^ can be considered as an e-perturbation of the symmetric 
matrix A. Thus the desired estimates can be obtained by the standard 
perturbation argument. 

Moreover, to prove the first inequality about the closest eigenvalue A 
one can use a more direct argument Q. Namely, for \i = A the inequality 
becomes trivial, while otherwise 

\\v\\<\\(A-^I)- 1 \\.\\(A-^I)v\\ = r \-m\. 

\fJL — A\ 



Thus 

lM-A|<M< E . 



Proof of Theorem 1.1. It suffices to show that for each k = 1,2, ... ,n a 



projection to the space of piecewise constant functions of the A;-th eigenvector 



of the linear operator L introduced in Lemma 1.2 satisfies the equality of type 
described in Lemma l=apr with e = O(l). Let P ra : C 1 ([0, 1] — ► R 1 , R n ) be 
a projection operator acting from the space of C 1 functions / : [0, 1] — ► R 1 
to R n defined as follows: 

(P n f)i:=W-l)/n). 

Then, according to Lemma |l.l| , the following equality holds true: 

A n P n Ek = n^k^nEk + £ 



for a vector £ £ R n such that ||£|| < O(l). Now the application of Lemma [L3 
finishes the proof. | 



2 Properly scaled left triangular matrices 

Results of the previous section might be generalized for a more general class 
of left triangular matrices. Let <f) : [0, 1] — > (0, 1] be a C 1 positive function. 
We define a family ofrexn left triangular matrices generated by the function 
cj) as follows: (Afyij '■= <P((j — i + l)/n) for i < j and zero otherwise. Clearly 
the binary left triangular matrices satisfy this property for <p = 1. 

Again similarly to the argument in the previous section one can prove the 
asymptotically linear dependence on n of the eigenvalues of the matrix A^ 
using the property that the main contribution to the "shape" of eigenvectors 
comes from the corresponding eigenfunctions of the integral operator L^ : 
f(x)^ti- x J>(s)f(s)ds. 

Theorem 2.1 Let jjl and /„ be respectively an eigenvalue and an eigenfunc- 
tion of the operator L^ (i.e. L^f^ = (xf^). Then for each positive integer n 
there exists an eigenvalue A = nfj, -\- 0(1) and an eigenvector e S R n of the 
matrix A^ such that (e)j = /«((i — I)/"-) + 0(l/n) for each i = 1, 2, . . . , n. 



Since in this more general case we are not able to give an explicit expres- 
sion for the asymptotic spectrum, we shall demonstrate also another (more 
direct) way to derive the asymptotic linearity of the spectrum. Assume that 
the function (j> is nonincreasing and denote by S 1 " the set of C 1 monotonically 
decreasing concave positive functions / : [0, 1] — ► [0, 1] with /(0) = 1 and by 
S™ C R™ - the set of vectors v such that v\ := f((i — l)/n) for some / £ Sf. 
We introduce also the following nonlinear operator B n : R n — > R ra U (oo)™ 
defined by the relation: (B n v)i := {Af l v)i/(Af l v)\. 

Lemma 2.1 The set S n is invariant with respect to the operator B n and for 
any vector dgS* the sequence {B™v} m converges as m —> oo in L 2 -norm 
to the leading eigenvector e{ of the matrix A$. 

Proof. Since the coordinates of the vector v are positive we get that for 
any i 

(A+v) i+1 - (Atv)i = -<f>((n - i + l)/n)« n _ i+ i < 0. 

On the other hand, due to the fact that the coordinates of the vector v 
decrease monotonically we obtain 

V) i+1 - 2(Aiv)i + (^v)i-l = <f>(- " )v n - i+2 ~ <f>(~ " )«n-i+l 

n n 

,n-i + l., . ( ,,n — i + 1. ,,n — i + 2. \ 

-){Vn-i+2-V n -i+l)-\ <f>{ ) ~ <f>{ )V n ~i+2 «n-i+2 



n \ n n 

fi — 2-1-1 

< <t>{ ){v n -i+2 ~ Vn-i+l) < 0. 

n 
These two inequalities prove the monotonicity and the concavity of Af^v 
respectively. The normalization by the positive number (Af l v)\ finishes the 
proof of the first statement of Lemma 2.1. Now noticing that the matrix A^ 
has only non negative entries we derive the second statement from the well 
known Perron theorem. ■ 



Thus choosing vectors v from the set S™ providing the smallest and the 
largest contribution to (A^v)i (namely (v)i = (n — i + l)/n and (v)i = 1 
respectively) we immediately obtain that 

n/2 / (1 - s)(f>(s) ds + 0(1) < A[ n) < n [ <f>(s) ds + 0(1). 
Jo Jo 

In a similar (but somewhat involved) way one can study also other eigen- 
values and eigenvectors of this matrix. 



3 Application to chaotic dynamics 

Now we apply above results for the analysis of ergodic properties of a certain 
family of chaotic maps. We refer the reader for the necessary definitions and 
statements for example to [||, |J. 

Fix a positive integer n and consider a piecewise linear map map T n : 
[0, 1] — ► [0, 1] from the unit interval into itself defined as 

Tn\[(i-l)/n,i/n)( x ) '■= ~% 

for all i £ {1,2,... ,n}. This map is topological^/ mixing and its second 
iterate is piecewise expanding (i.e. the modulus of the derivative of the map 
is larger than 1 for all points where it is well defined), therefore this map 
has a unique absolutely continuous invariant measure (see (l|). 

Lemma 3.1 The topological entropy of the map T n satisfies the equality 
hto P (T n )=ln(^ +0(1)). 

Proof. Notice that the map T n is Markov with respect to the partition to 
its intervals of monotonicity, i.e. it maps each interval of monotonicity to the 
union of some other intervals of monotonicity. Therefore one can construct 
the corresponding topological dynamics as a shift operator a n in the state 
of sequences y = {yi,y2, ■ ■ • } with the alphabet consisting of n symbols (i.e. 
j/i £ {1,2, ..., n} and (cr n y)j := yj+i), satisfying the property that for each 
k £ {1,2,... , n} the symbol k may be followed by any symbol from the set 
{1,2,... ,n — k + 1}, but not the symbol from the complement to this set. 
Thus the corresponding transition matrix A n [g] (matrix consisting of zeros 
and ones, describing possible transitions between symbols of the alphabet 
under the dynamics) is the left triangular binary nxn matrix. Now since the 
topological entropy ht op for both the Markov map T n and the shift operator 
o~ n is known to be equal to the logarithm of the largest eigenvalue of the 



transition matrix |2j] we derive from Theorem LI the desired relation for the 
topological entropy. | 
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